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ABSTRACT 

McLean  and  Saff man's  [1]  model  for  the  fingering  in  a  Hele-shaw  cell  is 
solved  numerically.  The  results  suggest  that  a  countably  infinite  number  of 
solutions  exist  for  non-cero  surface  tension.  This  set  of  solutions  contains 
the  solution  previously  obtained  by  McLean  and  Saff man  [ 1 ] . 
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SIGNIFICANCE  AND  EXPLANATION 


Two  dimensional  flow  in  a  porous  medium  may  be  studied  by  using  an 
analogue  proposed  by  Hele-Shaw.  The  analogue  is  based  on  the  fact  that  the 
mean  velocity  in  a  two  dimensional  porous  medium  and  the  velocity  of  the  flow 
between  two  parallel  plates  satisfy  the  same  equations* 

In  the  present  paper  we-  consider^the  steady  two  dimensional  flow  produced 
by  a  finger  advancing  between  the  two  plates.  The  analogous  porous  medium 

flow  occurs  in  oil  recovery.  This  problem  was  first  considered  by  Saffman  and 

i 

Taylor  [2] .  They  obtained  an  exact  solution  for  sero  surface  tension.  McLean 

and  Saffman  [1]  generalised  the  result  of  Saffman  and  Taylor  [2]  by  including 

the  effect  of  surface  tension  at  the  interface.  They  solved  the  problem 

numerically  and  obtained  one  family  of  solutions.  In  the  present  paper -we  ^ 

'  *  *  ! 
solve  the  problem, by  a  different  numerical  scheme.  Our  results  suggest  the 

existence  of  a  countable  number  of  solutions  for  non-sero  surface  teiwion. 

This  infinite  set  of  solutions  contains  the  solution  previously  obtained  by 

McLean  and  Saffman  (1]. 
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In  a  recent  paper  McLean  and  Saffman  [II  modelled  the  fingering  in  a 
Hele-Shaw  Cell  by  a  two  dimensional  potential  flow  with  surface  tension 
included  at  the  interface  (see  Figure  1).  They  introduced  the  notation  that 
U  is  the  velocity  of  the  finger,  2a  the  lateral  width  of  the  channel,  b 
the  transverse  thickness,  T  the  surface  tension  and  u  the  velocity.  In 
addition  they  denoted  the  ratio  of  the  width  of  the  finger  to  the  width  of  the 
channel  by  the  parameter  X. 

Taking  a  as  the  unit  length  and  ( 1  -  X )0  as  the  unit  velocity  they 
derived  the  following  nonlinear  integro-differential  equation  for  the  unknown 
shape  of  the  free  surface 

log  q(S)  “  -~I0  S^gTTg-  dS>  (1) 

KqS~(qs||)  -  q  «  -cos  0  (2) 

0(0)  -  0  q(0)  -  1  (3) 

0M)  »  "  \  q(  1 )  -  0  .  (4) 
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The  integral  in  (1)  is  of  Cauchy  principal  value  form,  the  variables  q  and 

A 

6  in  (5)  and  (6)  are  defined  in  terns  of  the  dimensionless  complex  velocity 

A  A 

u  *  iv  by  the  relation 

A 

-i0 


u  -  iv  *  qe  (8 

A 

the  flow  configuration  is  illustrated  in  Figure  1.  the  x-axis  is 
parallel  to  the  walls  of  the  channel  and  is  the  axis  of  symmetry  of  the 
finger.  The  points  A  and  B  in  Figure  1  correspond  to  s  »  0  and  s  *  1 
respectively. 

After  a  solution  for  6  and  q  is  obtained,  the  shape  of  the  finger  is 
given  by 


i6 


x(S)  +  iy(S)  "  -  ““  /g  dS  .  (9 


For  x  -  0,  Saffman  and  Taylor  [2]  obtained  the  following  closed  form 
solution 


lir-SIlMl 


(1-A)  +  S(2A-1) 


1/2 


do; 


8  “  cos  *q 


the  solution  (10)  leaves  the  parameter  A  undetermined. 

McLean  and  Saffman  [1]  diacretised  the  system  (1)  -  (4)  with  k  f  0  and 
solved  the  resultant  algebraic  equations  by  Newton's  method,  they  obtained 
one  family  of  solutions.  Romero  [3]  repeated  the  calculations  with  a 
different  numerical  scheme  and  different  initial  guesses.  He  obtained  two 
other  families  of  solutions. 

In  this  paper  we  present  a  systematic  numerical  procedure  to  compute  all 
the  possible  solutions  of  the  system  (1)  -  (4).  We  show  that  for  each  k  /  0 


there  exists  a  countably  infinite  number  of  solutions.  This  set  of  solutions 


contain  the  fan! lies  of  solutions  previously  obtained  by  McLean  and  Saffman 
[1]  and  Romero  [3].  Our  results  show  that  the  degeneracy  of  the  Saffman- 
Taylor  [2]  solutions  is  not  conp lately  removed  for  ic  /  0. 

Zn  order  to  find  all  the  solutions  of  (1)  -  (4),  we  define  a  modified 
problem  which  has  solutions  for  all  values  of  X  and  k.  This  modified 
problem  is  obtained  by  replacing  (4)  by 

q( 1 )  -  0  .  (11) 


Therefore  0(1)  becomes  a  free  parameter  to  be  found  as  part  of  the  solution. 

We  will  solve  the  modified  problem  defined  by  (1)  -  (3)  and  (11)  and 
obtain  solutions  for  all  values  of  X  and  k.  Then  we  will  obtain  the 
solutions  of  the  original  problem  by  selecting  among  the  solutions  of  the 
modified  problem  those  for  which  0(1)  “  -m/2. 

Following  McLean  and  Saffman  [1]  we  introduce  the  change  of  variable 

ST  -  1  -  •  (12) 


Here  T  is  the  smallest  root  of 


—  cot  wt  -  K 


(13) 


With  (12),  0  is  twice  differentiable  with  respect  to  c  at  both  end 

points.  McLean  and  Saffman  (1]  chose  y  ■  2  in  (12).  In  order  to  solve  the 
modified  problem  we  choose  y  m  4. 

To  solve  the  system  (1)  -  (3)  and  (11)  we  introduce  the  H  mesh  points 


-  (I-1)/N  I  -  1, . . . ,N 


(14) 


We  also  define  the  quantities 

0 


x  ■  ©  (J  1  “  ?j)  J  I  *  1,...,N  . 


(15) 


We  discretise  the  system  (1)  -  (3)  and  (11)  by  following  the  procedure 
outlined  by  McLean  and  Saffman  [1].  Thus  we  obtain  N  -  1  nonlinear 
algebraic  equations  for  the  N  -  1  unknowns  0  ,  Z  *  2,...,N.  For  given 


values  of  X  and  k  this  system  is  solved  by  Newton's  method 


In  Figure  2  we  present  numerical  values  of  6(1)  versus  X  for 
K  *  0.273.  As  X  approaches  zero*  0(1)  ♦  -*  and  the  finger  collapses  on 

A 

the  negative  x~axis.  As  X  approaches  one,  6(1)  oscillates  infinitely 
often  around  -t/2. 

Figure  2  shows  that  there  exists  a  countably  infinite  number  of  values  of 

X  for  which  0(1)  ■  -m/2.  The  solutions  corresponding  to  these  values  of  X 

are  the  solutions  of  the  original  problem. 

The  three  smallest  values  of  X  for  which  0(1)  ■  -  ~  are  0.54,  0.61 

and  0.67.  The  corresponding  profiles  are  shown  in  Figure  3.  The  value 

X  “  0.54  corresponds  to  the  solution  obtained  by  McLean  and  Saffman  [1]. 

Similar  results  were  found  for  other  values  of  k.  As  k  decreases,  the 

amplitudes  and  the  wavelengths  of  the  oscillations  in  Figure  2  decrease.  In 

Figure  4  we  present  the  values  of  three  smallest  values  of  X  as  functions  of 

„  1 

k.  The  three  curves  approach  X  «  ~  as  k  tends  to  zero.  The  lowest  of  the 
three  curves  is  the  family  of  solutions  previously  obtained  by  McLean  and 
Saffman  (1].  The  two  other  curves  correspond  to  the  families  of  solutions 
obtained  by  Romero  [3] .  The  values  of  X  presented  by  Romero  agree  with  our 
numerical  results  within  4%  for  k  >  0.2.  However,  Romero  did  not  obtain 
reliable  results  for  small  values  of  x. 
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